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ANALYTICAL  RELATIONS  BETWEEN  CONSTANTS  FOR 
GENERALIZED  VOIGT  AND  MAXWELL  VISCOELASTIC  MODELS 

ABSTRACT 

A  method  of  calculating  the  constants  of  a  Maxwell  model  from  the  con¬ 
stants  of  the  corresponding  Voigt  model  is  given.  The  derivation  is  based  on 
residue  theory  and  the  partial  fraction  expansion  of  the  transfer  functions. 
As  an  example,  constants  for  a  10-element  Maxwell  model  are  derived  from  the 
constants  of  a  10-element  Voigt  model  used  to  approximate  the  mechanical 
properties  of  polyisobutylene  at  25°C.  Formulas  for  the  creep  and  relaxation 
functions  are  also  given. 
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INTRODUCTION 


Voigt  and  Maxwell  models  are  frequently  used  in  analyzing  data  from 
mechanical  tests  of  viscoelastic  materials.  A  generalized  Voigt  model  is 
generally  more  convenient  for  representing  creep  and  complex  compliance  data, 
whereas  the  generalized  Maxwell  model  is  usually  used  for  relaxation  and  com¬ 
plex  modulus  data.  To  compare  the  results  of  different  types  of  tests,  a 
method  is  required  for  expressing  the  constants  of  a  given  type  of  model  in 
terms  of  constants  of  the  other  type.  Alfrey1  has  shown  that  relations  between 
models  consisting  of  three  or  four  elements  may  be  derived  by  elementary  alge¬ 
braic  methods.  When  the  models  consist  of  a  larger  number  of  elements,  the 
relations  among  the  constants  became  very  complicated,  and  a  method  based  on 
residue  theory  is  more  effective. 

On  several  occasions,  the  author  has  been  asked  to  discuss  this  problem 
and  outline  a  method  of  calculation  that  would  be  suitable  for  machine  calcu¬ 
lation.  For  this  reason,  a  detailed  analysis  of  the  problem  is  given  together 
with  a  numerical  example.  Gross  has  outlined  a  method  based  directly  on  the 
theory  of  limits.*  However,  his  argument  is  somewhat  difficult  to  follow.  It 
is  hoped  that  this  report  will  clarify  the  analytical  relations  between  Voigt 
and  Maxwell  models. 


ANALYSIS 

Springs  and  dashpots  are  the  fundamental  elements  of  a  viscoelastic  model. 
In  the  linear  theory  of  viscoelasticity,  it  is  assumed  that  the  stress  in  the 
springs  is  proportional  to  the  strain,  while  the  stress,  a  in  the  dashpots  is 
proportional  to  the  rate  of  strain,  The  stress-strain  law  for  the  entire 

model  may  then  be  expressed  as  a  differential  equation  with  constant  coeffi¬ 
cients  : 
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Viscoelastic  models  differing  in  configuration  are  considered  equivalent  if 
they  are  governed  by  the  same  stress- strain  law. 


^  Pages  6*4-65 ,  Reference  3 
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In  Eq.  (l),  the  low  order  terms  govern  the  behavior  at  long  times,  while 
the  high  order  terms  govern  the  behavior  at  short  times.  If  there  is  long¬ 
term  viscous  flow,  q^  =  0;  otherwise  >  0.  If  the  material  exhibits  instan¬ 
taneous  elastic  response,  xn  =  n;  if  the  elastic  response  is  retarded,  m  =  n  +  1. 
These  considerations  lead  to  the  four  classes  of  models  illustrated  in  Fig.  1. 
Voigt  and  Maxwell  representations  are  shown  for  each  class.  These  models  adhere 
to  the  topological  relations  given  by  Alfrey*  and  are  numbered  according  to 

Table  IV  of  Gross’s  monograph.**  In  this  figure,  the  E  represent  the  spring 

J 

constants  and  the  r\  represent  constants  for  the  dashpots . 

J 

The  limiting  properties  of  each  class  and  the  appropriate  stress- strain 
law  are  shown  in  Table  I.  It  is  clear  that  constants  for  a  given  model  can  be 
derived  from  constants  of  another  model  only  if  both  models  belong  to  the  same 
class  and  have  the  same  number  of  elements. 


Table  I.  Properties  of  Viscoelastic  Models  According  to  Class 


Long  Term 


Class 

Elastic  Response 

Viscous  Flow 

Form  of  Eq.  (l) 

I 

Instantaneous 

Absent 

m 

=  n; 

^0 

>  0 

II 

Instantaneous 

Present 

m 

=  n; 

«0 

=  0 

III 

Retarded 

Absent 

m 

*  n  + 

If 

o 

A 

c? 

IV 

Retarded 

Present 

m 

=  n  + 

If 

o 

(1 

c? 

*  Page  5^5 >  Reference  1 

**  Page  62,  Reference  J 


FIG.  i.  EQUIVALENT  VISCOELASTIC  MODELS 
CLASSIFIED  ACCORDING  TO  RESPONSE 
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Table  II.  General  Form  of  Viscoelastic  Functions** 


Viscoelastic  Function 

General  Form 

Complex  Compliance 

n+l 

^  ■  k  * 

Transfer  Function 

n+l  C, 

J(®)  -  Z  tA~  +  ji(b) 

J=3  8  1 

Creep  Function,  t  >  0 

n+l  \,t 

J(t)  -I  (1  -  e  J  )  +  e 
J=3 

(t) 

Complex  Modulus 

n  toi  *  E  • 

G*(io>)  =  V  — * — A  v  .  + 

G^la)) 

Transfer  Function 

n  C'  s 

0(6)  ‘  £kV=J +  °> 

) 

Relaxation  Function,  t  >  0 

n  l^t 

G(t)  =  I  C*  e  3  +  a  At) 
J-  K  3 

Table  III.  Remainder  Term  for  Each  Class  of  Model** 

Remainder  Terms, 

Class  of  Model 

Viscoelastic  Functions 

I 

ii  hi 

W 

J*1(  Ixo) 

1 

E1 

-i  +  — —  0 

Ex  in20) 
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1T)^0) 
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CD 
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A  +  -  0 

8 

B 

s 

€1(t) 

A 

A  +  Bt  0 

Bt 

E‘i 

0  E'^  +  iwri  2 

icon'2 

Gi(b) 

A' 

0  A*  +  B*s 

B»s 

A1 

0  A»  +  B'8(t) 

B*8(t) 

e  z  •  3  z 

**  Tn Tfebles  il  and  III,  <n  is  the  circular  frequency,  s  the  transform  variable 
in  the  Laplace  transform,  and  t  is  real  time. 
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The'  Voigt  and  Maxwell  representations  lead  to  expansions  in  terms  of 
partial  fractions.  Expressions  for  the  creep  and  relaxation  functions  are 
written  in  terms  of  the  constants  occurring  in  these  expansions.  The  general 
forms  of  the  various  viscoelastic  functions  are  shown  in  Table  II,  A  remainder 
term  depending  on  the  class  of  model  must  be  added  to  the  general  form  in  order 
to  obtain  the  complete  expression;  these  remainder  terms  are  given  in  Table  III, 


In  these  tables,  the  X,  are  roots  of  the  polynomial  equation 

u 


m 


qo +  qi8  +  •••  +  v5  =  0 


(2) 


while  the 


are  the  roots  of  the  equation 


p0  +  piS  +  ...  +  pnsn  =  0  (3) 

The  delta  function  6(t)  is  required  in  the  expression  for  the  relaxation  func¬ 
tion  for  Classes  III  and  IV,*  since  the  models  in  these  classes  cannot  yield 
instantaneously  to  a  load  of  finite  magnitude. 


Transformation  of  models  can  be  carried  out  only  for  models  belonging  to 
the  same  class.**  Consider  for  instance,  the  generalized  Voigt  model  and  the 
corresponding  Maxwell  model  of  Class  II,  and  suppose  the  constants  for  the 
Voigt  model  are  known.  The  roots  ...  of  the  equation 


A 


s 


=  0 


<*) 


are  calculated  by  some  method  of  successive  approximation.  These  roots  are 
real  and  negative,  and  are  interlaced  with  the  roots  X^,  X^  ... 


The  constants  C1 


j 


are  obtained  by  a  limiting  process. 


J(s) 


A  +  * 
s 


*  Page  ^/Reference  2 

**  Pages  6l  and  62,  Reference  3 
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bo  that 


But 


and 


C'  =  lim  (b  -  \i  ) 

G(s)  =-5Ti7 


c' 


J 


J 


The  indeterminate  form  on  the  right  is  evaluated  by  L* Hospital^  rule.*  The 
final  result  is 


CV 


~trW 

kk 


(5) 


A  + 


To  go  from  the  Maxwell  model  to  the  Voigt  model,  the  roots  X^,  X^  .  of 
the  equation 


La  - 


(6) 


are  calculated  in  the  manner  indicated  above.  The  coefficients  C.  are  calcu- 

J 

la  ted  from  the  formula 


C 


J 


l 


w; 

<ij  -  >*k)2 


*  This  rule  is  given  In  standard  texts  on  the  calculus . 
Page  J4 6,  Reference  4. 


(7) 


See  for  Instance, 
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The  constants  A  and  B  are  derived  from  the  relation 


G(s)  '  J(s)  =  1 


or 


=  1 


Let  s  0;  then 


B  = 


V  A.  ' 

L- 

The  value  of  A  is  found  by  allowing  s  to  become  infinite, 

1 


A  = 


* 


(8) 


(9) 


Formulas  for  the  other  three  classes  of  models  may  be  obtained  in  the  same 
manner. 


AN  EXAMPLE 

Consider  a  generalized  Voigt  model  used  to  approximate  the  dynamic  data 
for  polyisobutylene  at  25°C.  The  fitting  constants  were  calculated  by  a  semi- 
graphical  method  described  in  Reference  5,  The  calculated  values  of  the  complex 
compliance  generally  agree  with  the  measured  values  within  the  limits  of  experi¬ 
mental  error  over  the  entire  frequency  range. 
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COMPLEX  COMPLIANCE  OF  POLX ISOBUTYLENE 

25°C 

Measured  Values  (Ref  6)  Calculated  Values 

(Ref  5) 

Frequency 

J'(cu) 

J*'  (cju) 

J'(oj) 

cps 

cm^/ dynes 

cm^ /dynes 

cm  /dynes  cm  /dynes 

30 

204 

102 

205 

100 

40 

184 

93.6 

186 

100 

^5 

182 

107 

178 

100 

60 

156 

98.0 

159 

99.6 

80 

137 

97.8 

138 

97.7 

100 

119 

95-2 

122 

94.5 

140 

98.4 

88.1 

99.3 

86.7 

200 

79.0 

77-9 

79.3 

76.0 

280 

65.3 

65.0 

64.4 

66.2 

400 

51.2 

56.4 

51.4 

56.9 

600 

38.5 

46.9 

38.9 

47.3 

800 

32.1 

4l.O 

31-7 

40.8 

1000 

27.1 

36.4 

27.1 

36.3 

1400 

20.7 

30.1 

21.3 

30.3 

2000 

15.7 

26.5 

16.2 

25.O 

2800 

12.1 

20.5 

12.1 

20.5 

4000 

8.94 

16.6 

8.64 

16.1 

The  calculated 

values  tabulated  above  are  based 

on  the  following  model  constants 

MODEL  CONSTANTS:  10-ELEMENT  VOIGT  MODEL 

j  A 

B 

cj 

(cm^/dyne ) 

-9  2  / 

x  10  7  (cm  /dyn 

“9 

e  sec)  x  10 

2  -9 

(cm  /dyne  sec)  x  10  7 

rad/ sec 

3.365 

4.756 

1 

13110 

-165 

2 

71380 

-610 

3 

122600 

-3000 

4 

264500 

-15000 
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The  values  of  the  ^  for  the  corresponding  Maxwell  model  were  found  by 
J 

Newton’s  method.  The  fitting  constants  C*  were  then  obtained  from  Equation  5* 

J 

The  results  are  given  below. 

MODEL  CONSTANTS:  10-ELEMENT  MAXWELL  MODEL 


J 

1 

2 

5 

4 

5 


2  9 

(dynes/cm  )  x  10 

0.003498 

0.0009847 

0.004164 

0.01061 

O.2780 


“J 

rad/sec 

-17.52 

-222.5 

-1429 

-7611 

-151400 


The  roots  X.  and  ^  are  interlaced  along  the  negative  real  axis  in  the 
J  J 

manner  specified  by  the  general  theory;  moreover  the  constants  A,  B,  C.,  and 

J 


C 1 j  are  all  positive. 


The  values  of  the  individual  springs  and  dashpots 
derived  from  these  constants  will  be  positive ,  as  required  by  the  physical 
theory. 


The  relaxation  function  is  expressed  in  terms  of  the  constants  A,  B,  Cy 

and  \  for  the  generalized  Voigt  model.  The  creep  function,  on  the  other  hand, 
J 

is  expressed  in  terms  of  the  constants  C*  .  and  ^  of  the  generalized  Maxwell 

J  J 

model.  The  formulas  given  in  this  note  afford  a  simple  method  of  calculating 
the  mechanical  response  of  linear  viscoelastic  materials  under  a  variety  of 
loading  conditions  after  a  model  corresponding  to  a  particular  type  of  loading 
has  been  established. 


cl,  A, 


A.  S.  ELDER 
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